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Abstract 

A relativistic gas in a Schwarzschild metric is studied within the frame- 
work of a relativistic Boltzmann equation in the presence of gravitational 
fields, where Marle's model for the collision Operator of the Boltzmann 
equation is employed. The transport coefficients of bulk and shear viscosi- 
ties and thermal conductivity are determined from the Chapman-Enskog 
method. It is shown that the transport coefficients depend on the gravita- 
tional potential. Expressions for the transport coefficients in the presence 
of weak gravitational fields in the non-relativistic (low temperatures) and 
ultra-relativistic (high temperatures) limiting cases are given. Apart from 
the temperature gradient the heat flux has two relativistic terms. The first 
one, proposed by Eckart, is due to the inertia of energy and represents 
an isothermal heat flux when matter is accelerated. The other, suggested 
by Tolman, is proportional to the gravitational potential gradient and 
indicates that - in the absence of an acceleration field - a State of equi- 
librium of a relativistic gas in a gravitational field can be attained only if 
the temperature gradient is counterbalanced by a gravitational potential 
gradient. 



1 Introduction 

The research on relativistic gases by using the Boltzmann equation is an old 
subjcct in the litcrature. We can State that the Statistical description of a rela- 
tivistic gas began with the works of Jüttner in 1911 and 1926, when he succeeded 
to derive the equilibrium distribution functions for a relativistic gas that obeys 
the Maxwell-Boltzmann [TJ and the Fcrmi-Dirac and Bose-Einstein statistics 
[2]. The covariant formulation of the Boltzmann equation came later and was 
first proposed by Lichnerowicz and Marrot [3] in 1940. The establishment of 
the non-equilibrium distribution function and of the transport coefficients of a 
relativistic gas by using the Chapman-Enskog method was done by Israel [3] 
and Kelly [5] in the sixties of the last Century. These pioneers works were fol- 
lowed by several research papers in the literature concerning the study of gases 
in non-equilibrium states by using the Boltzmann equation in special rclativity. 
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However, there exist only few papers in the literature concerning the research 
of relativistic gases in the presence of gravitational fields within the framework 
of the Boltzmann cquation. The first works are due to Chernikov [üj [7] , who 
analyzcd the cquilibrium distribution funetions in some specific metric tensors 
and Bernstein [5], who determined the bulk viscosity for a relativistic gas in a 
Friedmann- Robertson- Walker metric. 

Some years ago the influence of gravity on the heat transport in a rarefied 
gas was analyzed in the works [ÜJ|TÖ] on the basis of a non- relativistic Boltzmann 
cquation, where it was shown that a gravitational field parallel (anti-parallel) to 
the temperature gradient increases (decreases) the heat transport. These results 
are related with the one derived by Tolman jlll IT2"] from a general relativity 
theory, who showed that a State of equilibrium of a relativistic gas in a gravita- 
tional field can be attained only if the temperature gradient is counterbalanced 
by a gravitational potential gradient. 

Reccntly a relativistic gas in a gravitational field by using the Boltzmann 
cquation was studied in the work |13j . where the dcpcndcncc of the heat flux 
on the gravity - the so-called Tolman's law - was obtained. In this work no 
metric tensor was assigned - although the components of the metric tensor 
appear explicitly in the cquilibrium distribution funetion - and the dependence 
on the gravitational field was connected with the Christoffel Symbol written in 
the Newtonian approximation. 

The aim of the present work is to analyze a relativistic gas in a Schwarzschild 
metric within the framework of the Boltzmann equation in the presence of grav- 
itational fields. The Marie model (Mj [15] of the collision Operator of the Boltz- 
mann equation is used and the transport coefficients of shear and bulk viscosities 
and thermal conduetivity are determined by the method of Chapman-Enskog. 
It is shown that the transport coefficients in the presence of a gravitational field 
are larger in comparison of their values in the absence of it. Furthermorc, the 
heat flux obtained has two relativistic terms, aside from the temperature gradi- 
ent. One of them - caused by the inertia of energy - represents an isothermal 
heat flux when matter is accelerated and was suggested by Eckart [16] . The 
other is proportional to the gravitational potential gradient and was proposed 
by Tolman [TTJ [T3] . It expresses that a State of equilibrium of a relativistic gas 
in a gravitational field and in the absence of an acceleration field can be attained 
only if the temperature gradient is counterbalanced by a gravitational potential 
gradient. 

The work is struetured as follows. In Section 2 the Marie model of the 
Boltzmann cquation, the equilibrium distribution funetion in the Schwarzschild 
metric are introduced as well as the definitions and balance equations for the 
particle four-flow and energy-momentum tensor. The non-equilibrium distribu- 
tion funetion, obtained by the use of Chapman-Enskog method, is the subject 
of Section 3. In this section the Euler equations of a relativistic gas in a gravita- 
tional field are obtained and the constitutive equations with the corresponding 
transport coefficients are determined for a viscous heat-condueting relativistic 
gas. In Section 4 we State the main conclusions of this work and in the appen- 
dices we introduce the isotropic Schwarzschild metric and show how to evaluate 
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the integrals when the distribution function depends on the components of the 
metric tensor. 



2 Boltzmann equation 

Let us consider a Single non-degenerate gas in a Riemannian space with line 
element ds 2 = g ßu dx^dx v where g ßv is the metric tensor. A particle with rest 
mass m is characterized by the space-time coordinates (x ß ) = (x° = ct,x) and 
by the momcntum four-vector = (p°, p), wherc c denotes the spccd of light. 

The constant length of the momentum four-vector g ßu p ß p u = m 2 c 2 implies 
that p° = Bi^mvL w here 

Po = \Jgoom 2 c 2 + {göigöj - goo9i ) p l p 3 ■ (1) 
The components of the four-velocity with U ß U ß = c 2 are 

tf" = (r c ,iv), r= 1 , (2) 

f \ ^ 

\ _|_ go» 
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and in a comoving frame v = yiclds U ß = ^y=j 0^ . 
For the isotropic Schwarzschild metric (|48|l we have 



P° = - Po = vS^VW, (3) 
5o 



where g = det(g pl/ ). 

In the phase space spanned by the space-time and momentum coordinates, 
the state of a relativistic gas is characterized by the one particle distribution 
function f(x a ,p a ) = /(x, p, t), such that f(x,p,t)d 3 xd 3 p gives the number of 
particles in the volume element d 3 x about x and with momenta in a ränge d 3 p 
about p at time t. 

The space-time evolution of the one-particle distribution function is governed 
by the Boltzmann equation. Here for simplicity we shall use a model equation 
for the Boltzmann equation, which rcplaces the collision term Q(f, f) of the 
Boltzmann equation by a collision model J{f). Two important model equations 
for relativistic gases were proposed by Marie [TU [15] and Anderson and Witting 
[TT] . The Marie model equation is given by 

T 

while the Anderson- Witting model equation reads 



J(f) = (/-/ (U} ), (5) 

T 

ig model equation reads 
J(f) = -^(f-f (0) )- (6) 



3 



In the above equations, r denotes a mean free time and /(°) the Maxwell- 
Jüttner distribution function. In the Marie model the Eckart decomposition 
for the particle four-flow and energy-momcntum tensor is used, whilc in the 
Anderson and Witting model the Landau-Lifshitz decomposition for thesc fields 
are employed. Note that in (|49p the four-velocity Ul refers to the Landau- 
Lifshitz decomposition. 

The collision terms of the model equations must fulfill the same properties 
as those of the true collision term of the Boltzmann cquation, namely 

1) For the summational invariant ip = p M , Q(f, /) and J(f) must satisfy the 
rclationship 

^Q(fJ)— =0, hence Uj{f) — = 0; (7) 
Po J Po 

2) The "H-theorem, or equivalently the tendency of the one-particle distribu- 
tion function to the cquilibrium distribution function, must hold 

/ QU, f) ^ /— < 0, hence / J(f) In ffl < 0. (8) 
J Po J Po 

For the proof of these properties one is referred e.g. to the book [T5j . 

Here we shall use the Marie model of the Boltzmann equation, which in the 
presence of gravitational fields reads 

('-*«>)■ 

Above, r* are the Christoffel Symbols and in a comoving frame the Maxwell- 
Jüttncr distribution function becomes 



thanks to ©2 and Qi. Here n,T and k represent the particle number density, 
temperature and Boltzmann's constant, respectively. Furthermore, ^(C) i s the 
modified Bossel function of second kind 

CV F(l/2) r.-tot.Z iN»-Va 



K "® = {i) wtmk «- ft (^-ir' äy, dl) 

which is a function of C = mc 2 /kT representing the ratio of the particle rest 
energy mc 2 and the thermal energy of the gas kT. The limiting case ( > 1 
corresponds to a gas in the non-relativistic regime at low temperatures, while 
(< 1 refers to a ultra-relativistic gas at high temperatures. 

The macroscopic State of a relativistic gas in a gravitational field may be 
described by the two first moments of the distribution function, the particle 
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four-flow N ß and the energy-momentum tensor T^ v . Their definitions in terms 
of the one-particle distribution function are given by 



d 3 p 
Po ' 



Po ' 



(12) 



The respective balance equations are obtained from the Boltzmann equation, 
yielding 



2V". 



0. 



T^, ß = 0, 



(13) 



where the semicolon denotes a covariant derivative. 

It is usual in the thcrmodynamic theory of rclativistic fiuids to decomposc the 
particle four-flow and the energy-momentum tensor in terms of quantities that 
appear in the theory of non-rclativistic fluid dynamics, namcly, particle number 
density n, energy per particle e, hydrostatic pressure p, non-equilibrium pressure 
w, pressure deviator ¥ ßV (the traceless part of the pressure tensor) and heat 
flux q ß . For Marle's model equation the decomposition used is that of Eckart 
(seee.g. HHCEB]) 



(14) 



1 P77 

T ßV = W -(p + w)A^ + — ( q »U u + q u U") + —U^U", (15) 



where is the projector 



(16) 



From the insertion of the Maxwell- Jüttner distribution function <jTÜ]) into the 
definition of the energy-momentum tensor (|12[) 9 and Integration of the rcsulting 
equation it follows that the energy per particle and the hydrostatic pressure 
read 

K 2 C. 

respectively. 



nkT, 



(17) 



3 Chapman-Enskog method 

In the Chapman-Enskog method the distribution function is written as / = 
f(°)(l + ip), where ip - the deviation from the Maxwell- Jüttner distribution 
function - is considered as a small quantity, i.e., \ip\ < 1. Furthcrmorc, the 
Maxwell- Jüttner distribution function is inserted on the left-hand side of the 
Boltzmann equation and the representation / = /(°' (1 + ip) on its right- 
hand side. By performing the derivatives it follows 



m 

T 



(/ - / (0) ) = / 



dn 
dx» 



V _ 
T 



1 - 



K 3 C , P T U T 



kT 



dT 
dx» 
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p^dUr c 2 dg l P y P \ _ \ __ Ar' t^f. 
' kT dx» 2kT dr U-p T 13 M r + kT 9l ij av U-p T 
c i o d 9a x x 1 \ m (0) 

-PP^r s ij—f^--f l fi- ( 18 ) 



2^kr r dr lJ r 

Hence, the deviation from the Maxwell- Jüttner distribution tp is determined as 
a function of gradients and derivatives of the metric tensor components. 

3.1 Euler's equations 

A relativistic fluid in the absence of the gradients of temperature and four- 
vclocity represents an Eulerian fluid. The determination of Euler's equations 
proceed as follows. First the multiplication of (fT5)l by \/—gd 3 p/po and Integra- 
tion of the resulting equation leads to the balance equation of particle numbcr 
density, namely, 

3 n 

UV d^ +nUV ' V=Q - (19) 
Next the multiplication of (|18[) by p^^—gd^p/po and subsequent Integration 
implies into an equation which is used to derive the balance equations for the 
energy density and momentum density of an Eulerian fluid. The cnergy density 
balance equation is obtained through the projection Up, yielding 

dT 

nc v U u — + pU»;„ = 0, (20) 

where c v is the heat capacity per particle at constant volume 
de , ( . 2 , .Ks K% 2 



- BT H C +5 ^ C "4 C V- (21) 

The momentum density balance equation results from the projection A^: 

mn^UW,, - ^p- -rnn^—^j-g^ = 0, (22) 
where we have introduced the gravitational potential 



GM . . <9$ GM . a?J 
$ = ---, wrth — = — Skj -. (23) 



We note that (|2"Tj) is a function of the ratio \Q(R)\/c 2 , which we can estimate 
at the surface of some bodies: 

1. Earth: M m « 5.97 x 10 24 kg; R m « 6.38 x 10 6 m; |$(i? e )|/c 2 « 7 x 1CT 10 ; 

2. Sun: M @ « 1.99 x 10 30 kg; i? Q « 6.96 x 10 8 m; |$(i? )|/c 2 « 2.2 x 1CT 6 ; 

3. White dwarf: M « 1.02M©; R m 5.4 x 10 6 m; |$(i?)|/c 2 « 2.8 x 1(T 4 ; 

4. Neutron star: M » M Q ; i? « 2 x 10 4 m; |$(i?)|/c 2 » 7.5 x 1CT 2 . 



6 



The above estimates imply that in most cases the approximation l'&l/c 2 <C 
1 is valid. Hence, the momentum density balance equation (|22[) in the non- 
rclativistic limiting casc ( > 1 and in the presence of a weak gravitational ficld 
reads 

-mn + (l + ^ + ...) ^=0- (24) 

Note that without the undcrlined terms (|24|) rcduces to the usual form of Ncw- 
ton's second law for a non-relativistic gas in the presence of a weak gravitational 
field. 



3.2 Viscous and heat-conducting relativistic gas 



In order to determine the energy-momentum tensor with the deviation of the 
Maxwell- Jüttner distribution function tp - given by (fT5|) - we first note that 
according to ([5]) and ([7]) the particle four-flow evaluated with the equilibrium 
and non-equilibrium distribution functions must have the same representation, 
i.e., 7V M = Ng = nU ß . From now on the index E will denote the value of a 
quantity at equilibrium, i.e., when evaluated with the Maxwell- Jüttner distri- 
bution function. However, in Marle's model the pressure and the energy per 
particle have different values at equilibrium and non-equilibrium states. This 
is not the case of the Anderson and Witting model equation §Q and of the füll 
Boltzmann equation (see e.g. [Tg]). 

In the following wc shall writc the constitutive cquations in a comoving framc 
where the componcnts of the projector become 



,00 



0, A' J >f ! 

.91 



1 



1 



1*1 
_v- 



(25) 



The knowledge of the energy-momentum tensor is obtained from the inser- 
tion of the representation of the distribution function / = + <p) - with 

<p given by (fT8)) - into its dchnition (fT2]) 9 and subsequent Integration of the 
resulting equation. 

Let us first analyze the projections 



ne = \ü ß U v T^ 



(26) 



which gives the following rclationships 



ro + (p- p E ) 



rkp 1 



3c^ (l 



1*1 

2c- 



Ke + Kl 



Kl 



7 



Kl 



or 2 3 

^2 



9CP 



(27) 



(e - e B ) 



Ä (l 



1*1 

2c- 



20 t4 + 3 Cß 



IS 



_ or 2 3 



017» 



(28) 



In the above equations it was introduced K E = -/^„(Ce)- 

We can relate the pressure difference (p — p^) with the energy difference 
(e — e^) by noting that pressure difference can be written as a temperature 



difference through (p — pe) = nk(T — Te) 



— 1/Ce)- If we expand 



the energy difference (e — e^) about Te and ncglcct terms up to second order 
in (T — Te) - since we are only interestcd in processes close to equilibrium with 
gradients of first order - it follows that e — ce — c E (T — Te). Hcnce these two 
last relationship together with (|26|) furnishes: 



P-Pe 



rp E k 2 



20^- 



3Ci 



13Ci 



Kl 
Kf 



2Ci 



9f 2 3 
JC E 



(29) 



Now we can obtain from (|27p and (|2"9"|) the final expression for the non- 
cquilibrium pressure 



w = —r\ 



dW_ 

dx3 ' 



(30) 



wherc 7/ denotes the coefficient of bulk viscosity, which is given by 



rk 2 p E 



3(c?) 2 (i-iay 



K E 
20^- 



1 #F 



+2 



cl 



cl 



2^, 2 
K. 



e^e + 3 Ce 




(31) 



The heat fiux and the traceless part of the pressure tensor are obtained from 
the projections 



A a A T A <Tr A 



T ß ", (32) 



and it it follows the constitutive equations 

dT T 



U a Ui, a + 



T 



9$ 



c 2 l-$ 2 /4c 4 0^ 



(33) 
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(6 ik ö jl + S a S jk ) 



-5 ij S k 



dx l 



(34) 



The coefficients A and /i are identified with the thermal conductivity and shear 
viscosity, respectively, and their expressions read 



A = 



ß = 



TB 



1 - 



2c 2 

t~Pe 



1 - 



1*1 

2c- 



k 



<e C + 5 



Kf 



K E ' 

^2 



El 



El 

K 2 E 



(35) 



(36) 



Equations (|30ß and (|34ß rcpresent the constitutive equations of a Newtonian 
relativistic fluid, while (f3"3")l the gener alized Fourier law. 

Let us analyze the transport coefficients given by (|31ß. (f35|) and (|36|) . All 
these coefBcicnts dcpcnd on gravitational potential through the component of 
the metric tensor gi(r) — (l + f^n:) and wc may infer from their expressions 
that the values become larger in presence of a gravitational field. The increase of 
the transport coefficients due to the gravitational field is rather small for stellar 
objects which are not too compact, i.e., for stellar objects where |*|/c 2 <C 1. 
In the absence of the gravitational potential the metric tensor reduces to the 
one of a Minkowski space-time. In this case gi = 1 and the expressions for the 
transport coefficients reducc to those of Marie jT5j HB] • 

In the Hmiting case of weak gravitational fields |*|/c 2 <C 1 and low temper- 
atures Ce ^> 1 the transport coefficients (f3"TT) . (|3"5j) and ([55)1 become 



V 



X 



= 5p g T 
6C 2 

ß = Pe t 

5kp E r 
~ 2m 



1 
1 

1 + 



21 
2& 

5 
2& 
3 



c 2 

m 

c 2 

21*1 



(37) 
(38) 
(39) 



The above expressions correspond to a non-relativistic gas in a weak gravita- 
tional field. 

The transport coefficients (|31|). (|35|) and (|36|) in the limiting case of weak 
gravitational fields |*|/c 2 <C 1 and high temperatures *C 1 read 



'/ : 

/' 

A 



54 
4p g r 

Ce 

^ 2 Pe> 
TCe 



1+1^ 

' 12 

r 2 

1 + ^ + . 



r 2 



?ln(- 
2 V 2 



41*1 



?7) 



21*1 



1 



21*1 



+ 



(40) 
(41) 
(42) 



9 



Here the expressions refer to a gas in a weak gravitational field and in thc 
ultra-relativistic regime. 

As a remark we call attention to the fact that in the literature of neutrons 
stars (see e.g. [ü]) the dependence of the transport cocmcients on the gravita- 
tional potcntial is not taken into account. 



3.3 Fourier's law 

Let us analyze Fourier's law with more details. It has three terms, namely, 

1 thc usual dependence of the heat flux in the gradient of temperature J^-; 

2 a rclativistic term proportional to —■^?U' J Uj; tT which rcpresents an isothermal 

heat flux when matter is acceleratcd. It is a small term that acts in a 
direction opposite to the acceleration and is due to the inertia of encrgy. It 
was proposcd by Eckart [16] within an irreversible thcrmodynamic theory; 

3 the third term ^ y 4e 4 -§§j is proportional to thc gravitational potential 

gradient. It was proposed by Tolman j!2) and indicates that in the ab- 
sence of the acceleration term a State of equilibrium of a relativistic gas in 
a gravitational ficlcl can be attained if the temperature gradient is coun- 
tcrbalanced by a gravitational potcntial gradient. In a weak gravitational 
field the equilibrium condition rcads 

~VT = -|, (43) 

which is the so-called Tolman's law with g being the gravitational field. 

By considering that the temperature field depends only on the radial coor- 
dinate T = T(r) and that there is no heat flux and acceleration field, it follows 
from ([33]) that the temperature field in a gravitational field must fulfill the 
differential cquation 



1 dT(r) 1 1 d$(r) 



T(r) dr c? 1 - $(r) 2 /4c 4 dr 



(44) 



The Solution of the above differential cquation for the boundary condition 
T{R) = 1 - where R is the radius of the spherical source - is 

T(r) = ) — — — y" , (45) 

(l + ^)(i-M^)' 

which is a decreasing funetion of r for r > R. In a weak gravitational field where 
|<I>(-R)|/c 2 is a small quantity (|45|) can be expressed as 

nr,-i-!T0-?)+W-?)'+-- <-) 
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It is worth to call attention that if we use the momentum density balance 
equation of an Eulerian fluid (|22[) in order to eliminate the acceleration term we 
get that the heat flux (|3"3")l becomes 



■ \ dT 



TK 2 dp 

nmc 2 K 3 dxi 



(47) 



We can transform the depcndcnce of the pressure gradient in (|47ß into a dc- 
pendence on the tempcrature gradient and particle numbcr density gradient 
through the use of the equation of State p = nkT and in this new description 
the presence of the gravitational potential gradient disappears. This indicates 
that the term related with the gravitational potential gradient in the heat flux 
is associated with the acceleration term. Howcvcr, in the work [13] the acceler- 
ation term was eliminated by the use of the momentum balance equation and 
Fourier's law has a dependence on the gravitational field and of the gradient of 
particle number density. We presume that the difference of the results is due 
to the fact that in the referred work there is no dependence of the equilibrium 
distribution funetion on the metric tensor components. 

4 Conclusions 

In this work we have analyzed a relativistic gas in a Schwarzschild metric within 
the framework of the Boltzmann equation in the presence of gravitational fields. 
The model of Marie for the collision Operator of the Boltzmann equation was 
used and the deviation from the equilibrium distribution funetion was deter- 
mined from the Chapman-Enskog method. Among other results it was shown: 
(i) apart from the temperature gradient the heat flux has two relativistic terms, 
one due to the inertia of the energy is connected with an acceleration, while the 
other is proportional to the gravitational potential gradient; (ii) the transport 
coefficients depend on the gravitational potential and become larger than their 
values in the absence of it; (iii) expressions for the transport coefficients in the 
limiting cases of low and high temperatures for weak gravitational fields were 
given; (iv) the radial temperature field of the gas is a decreasing funetion of the 
distance of the source of the gravitational field, when the heat flux and accel- 
eration field vanish; (v) the momentum balance equation for an Eulerian fluid 
in the presence of a weak gravitational field and in the non-relativistic regime 
reduces to the usual expression of Newton's second law. 

As a final comment it is worth to call attention that the present theory 
should not be valid in the case wherc $ — s- 2c 2 . In this case the the transport 
coefficients diverge and we are at the presence of the horizon of a Schwarzschild 
black hole. This limiting case will be the subject of a future investigation. 

A Appendix: Schwarzschild metric 

The Schwarzschild metric is the Solution of Einstein's field equation for a spher- 
ical symetrical non-rotating and uneharged source of the gravitational field with 
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total mass M. In terms of the spherical coordinates (r, 6, tp) it is given by 

^ = (l - 2 ^f) W - (i _Lm) ^ - + Sin2 ■ ( 48 ) 

where G is the gravitational constant. The isotropic Schwarzschild mctric (see 
e.g. [2Ü]) follows by introducing a new radial coordinate r 2 = öijX l x : > through 
the relationship r = r (l + §Mr) 2 so that (j48| becomes 

ds 2 = gr (r) (dx ) 2 - . 9l (r)Ö ij dx i dx j . (49) 
Above we have introduced the abbreviations 

(\ - GM) 2 f (ZM\ i 

^-J^Ef- 9 'M = ( 1 + S?)' (50) 

For the isotropic metric (|4"U)) the Cristoffel Symbols read 

roo = 0, r£. = o, o (i^j^k), r* 0j . = o, (5i) 

-.0 



~ 25! (r) dr ^ r ' ^ ~ 2 ff0 (r) dr * y r ' (52) 
whcrc the undcrlincd indices are not summed and 



(54) 



dg (r) _ 2GM (l - ) rf gl (r) _ 2gM / | GM\ 3 

dr c 2 r 2 (1 1 gm\ 3 5 dr c 2 r 2 \ 2c 2 r ) 

V 2c 2 r ) 

B Appendix: Evaluation of integrals 

To evaluate the following integral 

Z = V=9 [e-^ uX P^, (55) 
7 Po 

we choose a comoving frame so that the above integral reduces to 



Z = 4*J^$ r e -^/(^)|p| 2 ^, (56) 
Jo Po 

through the introduetion of spherical coordinates d 3 p = |p| 2 s'm8d\p\d6dip and 
by integrating in the angles < 9 < tt and < tp < 2ir. If we use the relationship 
^2 and change the variable of Integration through po = mey^/gö we have that 



l0 m 2 c 2 , o ,, mc ydy 

|p| 2 = (2/ 2 -l), d\p\ = — f » ~~T j (57) 



9i y/9i (y 2 - 1) 
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and (J5ÜJ) reduces to 

/oo 
e-^iy 2 - l)?dy = 4>irmkTK 1 (C), (58) 

i.e., Z is given in terms of the modificd Bossel function of second kind. 
For the evaluation of the integral 

Z 1 = y/=g\ — -, 59 

J U T p T Pq 

we change the variables and introducc 

mc 

Po = mcyfgöcosht, |p| = — =sinhi, (60) 

and gct 

Z x = Airm [°° cosh *- 1 e -Cco»ht dt = 4?rm _ Kh(0) , (61) 

Jq cosh t 

thanks to the dcfinition 

/oo roo — £coshi 

Ki n _ 1 (t)Ä = ^ cQgh „ - dt. (62) 

Following the same methodology it is casy to obtain 

Z^ = [ p » e -T^ UxpX — = ATtm 2 kTK2(QU», 

' J Po 



(63) 



Zr = ^/^^^=W^, (64) 

J U T p T po ( 

and to derive further exprcssions for Z^ 1 '" 1 *™ and Zf 1 
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